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Abstract 

The off-shell aspects of the one-pion-exchange potential (OPEP) are dis- 
cussed. Relativistic Hamiltonians containing relativistic kinetic energy, rela- 
tivistic OPEP with various off-shell behaviors and Argonne f is short-range pa- 
rameterization are used to study the deuteron properties. The OPEP off-shell 
behaviors depend on whether a pseudovector or pseudoscalar pion-nucleon 
coupling is used and are characterized by a parameter fi. We study potentials 
having //= — 1, and +1 and we find that they are nearly unitarily equivalent. 
We also find that a nonrelativistic Hamiltonian containing local potentials 
and nonrelativistic kinetic energy provides a good approximation to a Hamil- 
tonian containing a relativistic OPEP based on pseudovector pion-nucleon 
coupling and relativistic kinetic energy. 

PACS numbers: 21.30.-x, 24.10.Jv 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The modem high-quality two-nucleon potential models include three local potentials: 
Nijmegen II Reid93 and Argonne fig (AV18) 0, and two nonlocal potential: Ni- 
jmegen I and CD-Bonn . These potentials accurately fit the A^A^ scattering data of the 
Nijmegen database |^ and the deuteron binding energy, and are essentially phase-equivalent. 
The three local potentials predict very similar deuteron D-state probability Pd (5.70, 5.64 
and 5.76%) and triton energy (—7.63, —7.62 and —7.61 MeV) respectively, however the non- 
local CD-Bonn potential gives rather different results: 4.83% and —8.00 MeV. The Nijmegen 
I predicts similar Pd as the local Nijmegen II but slightly larger triton energy of —7.72 MeV. 
The experimental value of triton energy is —8.48 MeV, and Pd is not an observable. One 
may naturally ask: why is CD-Bonn so unique among the modern potentials? Is one model 
better than the other, or are all modern potentials equally good in predicting properties of 
nuclei after correctly including all other ingredients such as three-body forces, relativistic 
effects, meson-exchange currents, etc.? What role do nonlocalities play in properties of light 
nuclei and nuclear matter? 

To answer these questions, one has to look into the microscopic description of nu- 
clear interactions. Theoretically as well as experimentally it is known that the one-pion- 
exchange potential (OPEP) provides the long range part of the two-body potential v^^. The 
intermediate- and short-range parts of vnn are obtained by fitting the experimental data. 
The three local potentials have the same nonrelativistic on-shell OPEP at long range, but 
different phenomenological parameterizations at intermediate and short range. Nijmegen I 
has nonlocalities in the central part of the potential at short range. CD-Bonn, however, 
contains off-shell nonlocalities in OPEP predicted by pseudoscalar pion-nucleon coupling as 
well as nonlocalities in the phenomenological shorter range parts. The fact that the local 
potentials and the Nijmegen I model predict very similar deuteron properties which are very 
different from those predicted by CD-Bonn seems to suggest that differences in OPEP could 
be the main cause. 
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Inspired by the above comparison among the modern potentials, we want to understand 
various aspects of OPEP off-shell nonlocalities in this work. The nonlocality in the shorter- 
range part of the potential may also play some role, but it is not our primary interest here. 

There has always been some ambiguity about the off-shell effects in the two-body po- 
tentials. Friar |^ recently categorized these ambiguities into three types: (i) those caused 
by an energy- dependent potential which occur naturally when expanding energy denomi- 
nators in Schrodinger perturbation theory; (ii) those arising from unitary transformations 
of field variables used in the Lagrangians; (iii) those due to different choices of relativistic 
Hamiltonians. In this work we are only interested in understanding the second type of am- 
biguities. In Friar's notation, potentials having different off-shell forms are characterized 
by parameters fi and u, where the //-dependence comes from whether using pseudoscalar 
(PS) or pseudovector (PV) relativistic pion-nucleon interactions, while the i/-dependence 
comes from the retardation effects. Here we are only interested in the //-dependent off-shell 
behaviors, but neglect the retardation effects which seem to be relatively unimportant. In 
Friar's notation, we choose z/=l/2 (no retardation) throughout this work. 

The relativistic OPEP are identical on-shell, regardless of the assumed coupling (PS or 
PV), but they differ off-shell. There is no unique description for the off-shell behavior of 
relativistic OPEP. It depends upon a parameter /t whose common choices are /t=— 1 
(PS coupling) used in CD-Bonn |^, /i=0 (minimal nonlocality) used in our earlier work 0, 
and /t=+l (PV coupling) favored by conventional Chiral Perturbation Theory (CPT). As 
suggested by Friar [^|[ two decades ago, two-body potentials differing in the value of /i 
are related by unitary transformations up to order All these potentials are correct 

to this order, and predict identical observables, even though they have different forms. In 
this view, when various two-body potentials are combined with their corresponding three- 
body potentials, the triton energy should be accurately predicted; when electromagnetic 
current operators are treated consistently with the potential, the electromagnetic observ- 
ables should be independent of the choice of /t 0]. This is rather interesting and useful 
because if there indeed exists such unitary equivalence, then we could choose the simpler 
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and computationally-easier local potentials together with their three-body forces, and con- 
sistent electromagnetic charge and current operators to study properties of light nuclei. 

The purpose of this paper is to study deuteron properties with potentials containing 
yU=— 1, 0, +1 relativistic OPEP (Sec. II), and to examine the unitary equivalence of these 
potentials (Sec III). We also attempt to find local potentials which would give a smaller 
value for deuteron Pd, characteristic of the /i=— 1 nonlocal potentials, but with no success 
(Sec IV). The conclusions are given in Sec. V. Some of the detailed derivations involved in 
this work are given in an Appendix. 



II. OFF-SHELL EFFECTS IN DEUTERON 



Consider a relativistic Hamiltonian 



HR = 2^p^ + m'' ~2m + v^ + VR, (2.1) 
in the frame in which the deuteron is at rest. The transformation of this Hr to other frames 



is discussed in Ref. ||T0[. Here v-,^ is the OPEP containing off-shell nonlocalities, and vr is 
the remaining part of the potential which is phenomenological. In this work we use AV18 
parameterization for vr. As suggested by Friar after neglecting retardation effects, the 
relativistic OPEP can be expressed in the following general form 



<(P',P) 



f^NN 'Ti-Tj m m 



ml ml + 0^ E E' 
+fi X {E' - E) 



CTj ■ q CTj ■ q 



(Ti P (Tj p (Ti P (Tj P 



(2.2) 



E' + m E + m 

Here and m are the pion and nucleon mass respectively, f^NN is the pion-nucleon coupling 
constant, p and p' are the momenta of particle i in the center of mass frame before and after 
the interaction, q=p' — p is the momentum transfer, E=y/m''^ + p'^, and E'=y/m'^ + p''^. The 
/i-dependent term corresponds to nonlocalities in configuration space; it vanishes on-shell 
where E=E'. When used in momentum space to solve for deuteron properties, equation 
(172) is multiplied by the irNN form factor E{q) to ensure convergence 



In the present work we use the cutoff mass A^=1.2 GeV. 

The expression for yU=— 1, i.e., f^^~^ can be easily derived using second-order covariant 
perturbation theory with a pseudoscalar pion-nucleon couphng 

Hl%^ = iGij^^rMi (2.4) 

with G'^='irn?' fl^ / m^. When applying the same technique to the pseudovector-coupling 

HrNN = -^V^T'^Tsr.^a,^,, (2.5) 



we can not obtain the /i=+l expression The reason is that the derivative coupling 

of the pion field yields a term dependent on the pion energy go, and within second-order 
diagrams conservation of energy at each vertex requires qo=0. The expression corresponding 
to qo=0 has been given in Ref. |TI| up to order /m?, and is different from f^^^^ of that 



order. The correct expression for +^ [Eq. ( |2.2| )] can be obtained from fourth-order 
calculations |T2| . 



The /i=0 ("minimal nonlocality" ) corresponds to on-shell relativistic OPEP 

„M=0/ / _ flNN ^^ • q • • Tj mm 

mi mi + EL' 

TV TT ' ^ 

and has been studied in detail in Ref. 0. In the nonrelativistic limit, E ^ E' ^ m, and 
Eq. (|2.6| ) yields nonrelativistic local OPEP used in the three local potentials (Nijmegen II, 
Reid93 and AVI 8) 



m^ m^ + q^ 



When Fourier transformed to configuration space, v^^ [Eq. ( p^.7| )] becomes the familiar 
sum of the tensor and spin-spin interaction terms. The (5-function singularity present in the 
spin-spin term is usually removed or "smeared" by various form factors, e.g., exponential 
form factor in Nijmegen II and AV18, dipole form factor in Reid93 and CD-Bonn [Eq. ( |2.3| )]. 



We are interested in understanding the difference between the local and nonlocal potentials. 
It is sufficient to consider only one of the three available representations of a local potential, 
which we choose to be AV18. The results can probably be generalized to all modern local 
potentials. 

The OPEP in AV18 is given by the following expression in configuration space 

^.(r) = 1^.^ [YAr)<Ti ■ (Tj + T^{r)Sij] n ■ r,-, (2.8) 

where 

Sij — 3<Ti • rcTj r — (Ti ■ crj. (2-11) 

The rest of the AV18 is phenomenological and uses the Woods-Saxon and functions. 
When projected onto spin S=l and isospin T=0 channel, the overall potential is expressed 
as 

vs=i,T=o{r) = Voir) + Vt{r)Sij + vis{r)L ■ S + vi2{r)L^ + ^;,.2(r)(L • S)^ (2.12) 

where the five terms are called the central, tensor, spin-orbit, quadratic orbital angular 
momentum and quadratic spin-orbit terms, respectively. The OPEP is obviously included 
in the central term and is the main contribution to the tensor term. 

It is interesting to study deuteron properties using the isoscalar part of AV18 but replace 
its OPEP with v/^=-i'0'+i (for future reference, these potentials are denoted as Vi^~^'^'~^^). 
We construct the relativistic Hamiltonian containing Vig and relativistic kinetic energy to 
be phase equivalent to the nonrelativistic Hamiltonian containing the isoscalar part of AV18 
and nonrelativistic kinetic energy. The short-range part of the potential vr depends on /i 
implicitly in order to reproduce the phase shifts. We adjust the fifteen free parameters in vr 
and fit the phase shifts in S=l, T=0 channel, i.e., ^5*1, Ei, ^Di, ^D2, ^D^, E^, ^G^, and the 
deuteron binding energy with great accuracy. The deuteron properties are listed in Table I 
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and the deuteron wave functions are plotted in Fig. |l[ The first row of Table I shows the 
experimental results; it is followed by results obtained with five modern potentials, and the 
last three rows are the results obtained with f fo"^"*^'"'"^^. 



state probability (Pd=5.73%) to that of the nonrelativistic AV18 (5.76%). Their deuteron 
wave functions are also similar as shown in Fig. |l]. The f ig"^ binds the triton slightly more 
(~0.1 MeV) than AV18 [§. These similarities between fi8~° and AV18 seem reasonable 
because the Hamiltonian is constrained such that the nonlocalities from the f i8~° and rela- 
tivistic kinetic energy Tr have to cancel exactly in the deuteron to give the same binding 
energy. The small variations in and triton binding energy are due to the almost perfect 
cancellation between f and Tr nonlocalities. 

In the cases of /i=— 1 and +1, there does not seem to have such a perfect cancellation. 
The deuteron Pd predicted by ffs""^ is 4.98%, noticeably smaller than that of ffg"^ (5.73%), 
while the of fig""^^ is 6.26%, noticeably larger than that of v^^^ . This is entirely due to 
the off-shell terms in and absent in because f j^g""^'"'"*"^ all have the same 

short-range parameterization of AV18; only the values of parameters are slightly different 
in order to fit the data. The Pr, of f i8~~^ (4.98%) is very close to that of CD-Bonn (4.83%) 
as expected since both potentials have the same OPEP "w^^^^ apart from different pion 
cut-off mass used in the form factor [Eq. ( p.3|) ]. We use 1.2 GeV for A^ while CD-Bonn 
uses 1.7 GeV. When we use A^=1.7 GeV in our v^^~^ as a test and refit the phase shifts, 
the deuteron properties remain almost the same (Pd=4.97%, D / S ratio=0.0240, (5d=0.0267 
fm^). Therefore the differences between ffg""^ and CD-Bonn must come from those in the 
short-range parts. This is further confirmed by the similar D-state wave function at long 
range (r > 1.5 fm). The difference in the deuteron quadrupole moment given by the two 
potentials is ~1.5%. It is smaller than the 3% difference in P^. This can be understood 
from the following arguments. In the impulse approximation, deuteron Qd is given by 




has been discussed in detail in Ref. 0. It produces very similar deuteron D 




(2.13) 
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and it is an "outside" quantity whose major contributions come from the longer-range part 
of the wave functions which is fairly model- independent. On the other hand, the D-state 
probability Pu defined as 



oo 

2 



= / w\r)dr, (2.14) 



is an "inside" quantity lacking the extra factor in the integrand. Therefore is sensitive 
to the interior region of the wave functions where larger differences exist among various 
potential models. 

Figure also shows another interesting point. The AV18 wave functions are in-between 
those of f i8~^^ and v^^^^ . The f i8~~^ wave functions are farther away, and those of CD-Bonn 
are the farthest. The difference between the wave functions of AVIS and fig"^^ favored by 
CPT is rather small. This can be easily understood from the following first-order perturba- 
tive estimates. The largest contributions of OPEP come from coupling states with small p 
to large p' . Therefore we consider the extreme case where p=0, p'=q, and expand Eq. ( p.2|) 
in powers of (f' jvn? . This gives 

.r-(q,o)=.r(i-^ + ---), (2.15) 

<^°(q,0)=.r(l-^ + ---), (2.16) 

t'r-^^(q,0)=.r(l-^ + ---). (2.17) 

If we treat OPEP as a perturbative term in the Hamiltonian, the first order correction to 
the amplitude of a state with large momentum q is given by 

Here can be f or depending on whether the Hamiltonian is nonrelativistic or rel- 
ativistic. In the nonrelativistic case, E[(i) — E{0) ^ q'^/m, while in the relativistic case, 
E{q) — -E'(O) ~ g^/m x (1 — /Am?). Therefore the difference between the first order, large 
q, amplitudes of and is given by 
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r=-\ci) - ^^^(q) ^ -^^NR(q) (2.19) 

r=\cO - V^^^(q) ^ -^^^^(q) (2.20) 
^'^=+i(q)-^^^(q) f^iO. (2.21) 

This estimate suggests that the vl^^~^^ wave function is the closest to the nonrelativistic 
wave function, that of f^^'' is the second closest, and that of is the farthest. In re- 

ality, other components of v{p', p) contribute as well, and the OPEP can not be treated 
perturbatively. The exact calculations show (Fig. P that the nonrelativistic wave func- 
tion is in-between those of fig"^^ and fi8~°. This implies that the wave functions obtained 
from nonrelativistic local potentials such as AV18 and nonrelativistic kinetic energy provide 
excellent approximation to those with f ig"^^ and f and relativistic kinetic energy. 

The expectation values of the kinetic energy T, OPEP v^^ and the phenomenological 
short-range part vr of various potentials are compared in Table II. As we can see, the indi- 
vidual terms change slightly among various potentials, but the sum of them gives identical 
total energy of —2.242 MeV by construction. The experimental value of deuteron bind- 
ing energy (—2.224 MeV) is obtained from the full AV18 with electromagnetic interactions. 
From Table II, it appears that the individual terms of f fg"^^ are the closest to those of AV18, 
those of fi8~° are the second closest, and those of ffs""^ are the farthest. Also listed in this 
table are the on-shell (yo^i-shcu^ off-shell (^off-sheii^ contributions to the OPEP. The sum 
of the two terms gives the total [Eq. ( p.2|) ]. For AVIS and f i8~°, the off-shell contributions 
are obviously zero, and the on-shell contributions are the nonrelativistic OPEP [Eq. ( p.SI) ] 
and on-shell relativistic OPEP [Eq. ( |2.6| )], respectively. For v^g^^, the off-shell contribution 
is much smaller than the on-shell one because the former is of order /m? of the latter in 
the leading order. 

The deuteron charge form factors F"^ and T20 calculated in impulse approximation are 
plotted in Figs. and ^. Notice that the results predicted by AV18 are in-between those of 
f fg"^^ and f 1^°. The difference between the T20 obtained with CD-Bonn and AVI 8 is much 
larger than that calculated with f fg of various choices of /i implying that the nonlocality in 
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the shorter-range part of CD-Bonn is the main cause. 



III. UNITARY EQUIVALENCE 

Friar |^,^ suggested two decades ago that the physics of pseudovector and pseudoscalar 
couphngs of pions to nucleons is related via the equivalence theorem, originating from the 
Dyson transformation |]T3. Basically, the theorem states [Bl that the one-pion-exchange 



currents, together with matrix elements of the impulse charge operator calculated with a 
potential including OPEP, must be the same for both PS and PV couplings to order G^. 
As an interesting byproduct of Friar's work, the unitary transformation involved can change 
the deuteron D-state probability, indicating that this quantity is not measurable fl^ . 

It is interesting to test this theorem using our relativistic Hamiltonian [Eq. for 
/i=— 1, 0, +1. The Hamiltonian also contains a short-range potential vr which is implicitly 
dependent on /i in order to fit the data. We can rewrite Eq. (|2.1|) as the following 

HL^t = Tr + + A< + < (3.1) 

where Tr=2a/p^ + — 2m is the relativistic kinetic energy, is the on-shell relativistic 
OPEP given by Eq. ( p^.6D and Af^ is the off-shell term in Eq. ( |2.2| ) given by 



At;^(p,p)--/.x— ^-^--(E -E) [-^r^ eT^J- ^^-^^ 

In light nuclei, the small binding energy comes from the large cancellation between the 
kinetic energy and two-body potential. Hence and in Eq. ( |3.1| ) are of the same 



order of magnitude as Tr. The off-shell term Af^ is much smaller, its leading term is of 
order /m? of It seems that in deuteron the off-shell term Af^r is comparable to vr 

as shown in Table II. This is because deuteron is a loosely bound state with a very large 
rms radius, thus the OPEP accounts for more than 95% of the two-body potential, yielding 
a very small vr. In light nuclei with A > 2, the OPEP still accounts for a large portion of 
the two-body potential (>70%), but vr is about the same order of magnitude as the kinetic 
energy and OPEP |jl5| and At;^ is probably much smaller than vr. 

10 



The Hamiltonian [Eq. ( |3.1| )] can also be expressed, up to order /m?, as the following 
unitary transformation 

iJ^^V = e-'^i/oe^'' ^Ho + i [H,, U] , (3.3) 
where Hq can be conveniently chosen as 

Ho = i^eTact =Tr + + V^'. (3.4) 

The unitary operator which satisfies Eq. ( |3.3|) is found to be 

y'.NN r -r rn rn / a, ■ p- a, ■ p- _ ' P ' p\ _ ^3 5) 

^ 2 ml ml + q^ E E' \ E' + m E + m J ^ ^ 

It can be easily verified that the commutator 

[T^,^f/] = A<. (3.6) 
The wave function is transformed, consistently with the Hamiltonian, in the following way 

\r^' = e-'''\^o)^il-^U)\^o) (3.7) 

where {ipo) is the wave function of Hq obtained by solving Schrodinger equation. The meth- 
ods used to solve Eq. (|3^) for the deuteron are given in the Appendix. 

Earlier works by Friar Desplanques and Amghar and Adam et al. used a 
nonrelativistic kinetic energy /m in the Hamiltonian, and kept the leading order {jp' /rn^) 
terms in the "relativistic" OPEP and the unitary operator. This seems to be inconsistent 
with the relativistic considerations, nor is it computationally easier than the exact treatment 
since all the calculations involving the OPEP off-shell term are performed in momentum 
space. In this work we use a relativistic kinetic energy and keep the relativistic OPEP and 
unitary operator iU in their exact forms. 

From Eqs. ( |3.3| ), ( |3.4| ) and ( |3.5| ), we find that the unitarily transformed Hamiltonian is 



not exactly the same as the if^ 



exact 



Hl, = Hl^, + /\H^. (3.8) 
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The correction term is 



The two commutators [v^ iU] and 



<=°, iU 



+ 



- <• (3.9) 



are of the same order as Aff^, but the first 



commutator has the range of two-pion exchange interaction while the second one has even 
shorter range. Both can be absorbed into the phenomenological short-range part to 
obtain v^. If this is exactly the same as f^, then H^^~H^^^^^, and is unitarily equiv- 
alent to However, it is unlikely that v^=v^ as demonstrated in the wave functions 
shown in Fig. |^ and deuteron properties listed in Table III. 

As we can see in Fig. ^, the unitarily transformed deuteron wave functions are slightly 
different from the exact wave functions obtained by solving the Schrodinger equation, mostly 
at short range. The unitary transformation works rather well for the D-state. It also works 
well for the S'-state at long range but somewhat breaks down at very short range (r < 0.3 fm). 
This can be understood from the following arguments. The unitary transformation discussed 
in this work is related to OPEP. The long-range parts of the wave functions are mostly 
determined by the OPEP, thus they are related very well by the unitary transformation. 
The D-state is most sensitive to the tensor force whose main contribution comes from the 
OPEP. The S'-state, however, is sensitive to both central and tensor parts of the potential at 
short range, thus the unitary transformation starts to deviate from the exact S-state wave 
function at r < 0.3 fm. This deviation at such short range may have very small impact on 
the observables. For example, the deuteron binding energies calculated with the unitarily 
transformed wave functions for /i=— 1 and +1 are reduced by only ~8% as compared with 
the exact result, indicating that the unitary transformation is a good approximation to the 
exact solution. 

The deuteron properties obtained with unitary transformation are compared with those 
from solving Schrodinger equation in Table III. The differences between the exact solutions 
and the results of unitary transformation are obviously caused by the short-range part of 
the potential in A/7. These differences are rather small indicating that the potentials ufg 
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for /i=— 1, and +1 are, to large extent, unitarily equivalent. Desplanques and Amghar 



studied unitary equivalence between Paris (yU = 0) and the older Bonn (/i=— 1) potentials 
before the birth of the modern potentials. The differences between the results of Bonn and 
the unitarily transformed Paris, or Paris and unitarily transformed Bonn seem to be larger 
than the present work, indicating that the Paris and Bonn potentials are not as unitarily 
equivalent as our f models. The reason is probably due to the rather different structures in 
the short-range part of the two potentials, and possibly due to using the leading order terms 
in kinetic energy, OPEP and unitary operator instead of their exact expressions. Adam et 
al. [|ll|] did a similar study as the present work, except that they used a pure OPEP model 
to describe the two-nucleon system, and they also used leading order expressions for kinetic 
energy, OPEP and unitary operator. They studied deuteron properties for various values of 
fi and u. Their deuteron and Qd for z/=l/2 (no retardation) and /i=— 1, 0, +1 seem to 
be larger than ours, probably because they used a pure OPEP model instead of a realistic 
model which fits A^A^ phase shifts as well as deuteron binding energy. However, the amount 
of changes from before and after unitary transformation are similar to ours. 

The deuteron of current work is plotted in Fig. ^. The P^ obtained with unitary 
transformation seems to be linear in fi, while the exact results are slightly nonlinear. This 
can be easily understood from the expression of the unitarily transformed wave functions 
given in Eq. (|A16|) . We can rewrite the wave functions as 

RT (p) = (p) + /i X 6R2 (p) , (3.10) 

where Ro{p) and R2{p) are the normalized S- and -D-state radial wave functions, 6Rq{p) and 
(5i?2(p) are the second terms in Eq. ( |A16|) . In the dominant region of the wave functions, 
i.e., p < 5 fm~^, Ro{p) ^ 6Ro{p), R2{p) ^ 6R2{p)- The deuteron D-state probability is 
given by 

r roc roo 

p-" = pO -2P° / Roip)6Roip)p^dp + 2{l-P^) R2{p)5R2{p)p'dp 
L "'0 Jo . 

-{- terms of order /i^ and higher, (3-11) 
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where is the deuteron D-state probabihty of /i=0. The terms of order /i^ and higher are 
neghgibly small and can be neglected, therefore Pn obtained from the unitary transformation 
is approximately linear to fi. The exact Pd is slightly nonlinear reflecting the difference of 
the short-range part of the potentials. 



IV. LOCAL POTENTIALS 

Before we draw any conclusions, let's ask a question: is the small Pd a unique feature 
from or the large Po a unique feature from f ^^"'"^ ? In other words, can we find a 

local potential which fits experimental data and still gives Pd < 5% or Pd > 6%? 

To explore such a possibility, we attempt to lower the Pd to ~5% by using AV18 as trial 
local potential in our Hamiltonian, together with a nonrelativistic kinetic energy. There 



are fifteen free parameters in the AV18 representation [Eq. ( |2.12| )]. Just by varying these 
parameters to fit the data and simultaneously constraining Pd ~ 5% can result in good fits 
to ^5*1, ^Di, ^D2, ^D^, and ^G-^ in addition to deuteron binding energy, but not to the 
mixing angle Ei. At higher energies (200-400 MeV) Ei comes out significantly below the 
Nijmegen values. 

Since Pd is sensitive to the tensor force, in order to get smaller Pd, we need to reduce 
the tensor force. We have tried various approaches summarized below: (i) Multiply the 
entire tensor force Vt{r) by a factor (1 — e~'^'^^^), and adjust the parameter ci in addition 
to the fifteen parameters; (ii) Add a factor C2r(l — e~^'^^') to the tensor force, and vary 
the parameters C2 and C3 in addition to the fifteen parameters; (iii) Use the wave functions 
obtained with Vig^^ whose Pd is 4.98% and solve the Schrodinger equation for local Vc 
and Vt (denote them as "local /i=— 1"), keeping vis, f/2 and Vis2 fixed. These potentials are 
compared with those in the AV18 in Fig. |^. As we can see, the "local /i=— 1" Vc and Vt 
are both reduced at the short-range to reproduce the deuteron binding energy and obtain a 
smaller Pd- We then multiply or add various functions to Vc and Vt to reproduce this effect. 
None of these approaches yield a good fit- The largest discrepancies are for Ei. 
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This is not surprising because Ei is most sensitive to the tensor force, unhke any other 
states. When restricting to be ~5%, we are reducing the tensor force which in turn 
reduces the Ei. Other states can make trade-offs among Vc, vis, ■ ■ ■■ In contrast, when 
the OPEP tensor force is modified by the off-shell term in Eq. ( p.2| ), only the off-diagonal 
elements v{p,p') is reduced, but the Born term v{p,p) remains the same. This seems to 
suggest two things: (i) the tensor Born term plays an important role in fitting to Ei phase 
shift; (ii) weakening the tensor force off-diagonally does not seem to affect Ei, however it 
produces a smaller Pd of the deuteron. 

In summary, we couldn't find a local potential which gives a deuteron Pd as small as 5%. 
We can probably extrapolate this statement to Po as large as 6%. Therefore, it seems that 
the different deuteron Pd's are associated with the various off-shell behaviors of the OPEP. 

V. CONCLUSIONS 

We studied deuteron properties with potentials containing /i=— 1, 0, and +1 OPEP off- 
shell behaviors. We tested the unitary equivalence among these potentials by comparing to 
the exact solutions of Schrodinger equation, and we found the following interesting results: 

(i) The off-shell term involving the coupling of nucleon spin and momentum is the primary 
cause for the smaller deuteron P^ for CD-Bonn and ffg""^, and consequently a higher triton 
binding energy. We also found that it is not possible to get a deuteron with Pd as small as 
5% with the local OPEP [Eq. (U)]. 

(ii) The deuteron wave functions, charge form factors and T20 for /i=+l favored by CPT 
are very close to those of the nonrelativistic AV18. 

(iii) The t>fg"~^'°'^^ models are, to large extent, unitarily equivalent. 

(iv) The differences due to unitary transformations related to OPEP are smaller than 
the difference between AV18 and CD-Bonn potential. 

The charge form factors and T20 shown in this work are calculated in impulse approxi- 
mation. After including the meson-exchange current operator obtained with unitary trans- 
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formation consistent with the Hamiltonians and wave functions, the deuteron A and B 

turn out to be rather similar for f fg" 



structure functions as well as Ton turn out to be rather similar for f and AV18, but 



they differ significantly from those of CD-Bonn ||T^. This suggests that AV18 and CD-Bonn 
are not well-related by the kind of unitary transformation discussed in this work, i.e., the 
one that deals with OPEP. 

In summary, we find that a nonrelativistic Hamiltonian containing local potentials (AV18, 
Nijmegen II, Reid93) and nonrelativistic kinetic energy provides an excellent approximation 
to a relativistic Hamiltonian containing a relativistic OPEP of pseudovector pion-nucleon 
coupling and relativistic kinetic energy, in predicting properties of nuclei. This may matter 
little for two-nucleon systems such as the deuteron discussed in this work, or three-nucleon 
systems such as the triton discussed in Refs. 0,0 for which momentum space as well as 
configuration space computational methods can be applied. However, when we go up to 
nuclei having A > 3, momentum-space techniques face great computational difficulties, 
leaving the configuration-space calculations as the main approach for high-accuracy compu- 
tations P^JI^I, in which case a local nonrelativistic Hamiltonian is to be preferred. 
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APPENDIX A: UNITARY TRANSFORMATION OF THE DEUTERON WAVE 

FUNCTIONS 

Consider the case of deuteron. The Hamiltonian is given by Eq. ( |3.4| ) for /i=0. The 
wave function in momentum space is expressed as 

iV^o(p)) = i?o(p)3^iio(p) + R2{p)yii2{p) (Al) 

where Ro{p) and R2{p) are the S- and D-state radial wave functions, 3^l5j(p) are the spin- 
angle functions 

yLSJ = (-0"^ E(^. Ml = Mj - Ms, S, Ms\J, Mj)YL,Mr,\S, Ms) (A2) 

Ms 

Here J, S, L are the total, spin, and orbital angular momentum, respectively, and 
Mj, Ms, Ml are their projections along ^-axis. From Eq. (|3.7| ), the wave functions of 
/i=±l are given by 

mpr=^'r' = \mp)) - J ^,^u{p,p')\Mp')) m 

where the unitary operator is given by Eq. (p.5|) multiplied by the form factor [Eq. (|2.3| )]. 
We can rewrite 

ml + q^[Al + q^)~ml + q^ + + ^ - "^^^ dAl Al + q^' ^ ^ 

Then use 

-1^2 =^T.Qli^mjYLimn.{p) (A5) 

mi + q^ pp' 

where Qi{zm^) is the Legendre function of the second kind and 

_ + p'^ + ml 

Because deuteron has only /=0 and 2 states, we will only need Qq and Q2 

Qq[z) = -In 



(A6) 



2 \z-l 

Q2{z) = ^{3z' - l)Qo{z) - h (A7) 
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Using these, Eq. ( |A4| ) becomes 



where 



PP Im 



d 



Qi{p',p) = QiizmJ - Qi{zaJ + {Al - mD^QiizA^^ 



Then we use the following relations 



(A8) 



(A9) 



CTi ■ p era ■ p = -p^ (5'i2(p) + CTi ■ 0-2) 

^12(P)3^110(P) = V^3^112(P), 
512(P)3^112(P) = V83^110(P) - 23^112 (P) 

CTi ■ cr23^iio(p) = 3^iio(p), 

CTl ■ Cr23^112(p) = 3^112 (P)- 



(AlO) 
(All) 
(A12) 

(A13) 
(A14) 



Finally, putting everything together, we get 



i^'^(p))-' = Rrip)yuo{p) + RTWiuiP) 



(A15) 



where 



Rrip) = Roip) + /i X / ^ {{E' - E)Qo{p',p)Ro{p') 



i?r(p) = ^2(p)+/ix 
+ 



Svr^ ml EpJ E' 
{E' - m)Qo{p',p) -{E- m)Q2{p\p) 

1 fnNN^^ fP'dP 



R, 



^2 ^2 J J^i 

{E' - m)Qo{p',p) -{E - m)Q2ip',p) 



8{E'-E)Qo{p',p)Ro{p') 



R, 



(A16) 



The wave functions should be normalized to 1, then Fourier transformed to configuration 
space to obtain the S- and D-state wave functions shown in Fig. HI 
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FIGURES 

FIG. 1. Deuteron S- and D-state wave functions. 

FIG. 2. Deuteron charge form factors squared for spin projections M(i=0 and 1 in impulse 
approximation, for various potentials. 

FIG. 3. Deuteron tensor polarizing power T20 in impulse approximation, for various potentials. 

FIG. 4. Comparison of deuteron wave functions obtained from unitary transformation and 
those from solving Schrodinger equation. 

FIG. 5. Deuteron D-sate probability as a function of fi. The circles represent the Pd ob- 
tained by solving Schrodinger equation, while the diamond symbols represent those from unitary 
transformation. 



FIG. 6. Comparison of the "local /x=— 1" central and tensor potentials with AV18. 



21 



TABLES 



TABLE I. Deuteron properties predicted by the modern potentials and wfg for fJ,=—l, 0, +1. 





Character 


Qd (fm^) 


D/S ratio 


Pd (%) 


Experiment 


nonlocal 


0.2859(3) 


0.0256(4) 




AV18 

Nijmegen II 
Reid93 
Nijmegen I 
CD-Bonn {n=-l) 


local 
local 
local 
nonlocal 
nonlocal 


0.270 
0.271 
0.270 
0.272 
0.270 


0.0250 
0.0252 
0.0251 
0.0253 
0.0255 


5.76 
5.64 
5.70 
5.66 
4.83 


M=-l 
^18 

u=0 
^18 

^18 


nonlocal 
nonlocal 
nonlocal 


0.266 
0.271 
0.272 


0.0253 
0.0258 
0.0260 


4.98 
5.73 
6.26 


TABLE II. Expectation values of kinetic energy T, OPEP 
short-range part vr of various potentials (in MeV). 


and the phenomenological 




AV18 


/i=-i 


^i=o 

Vl8 


M=+l 


{E}^ 


-2.242 


-2.242 


-2.242 


-2.242 


(T) 


19.882 

-21.355 
-0.769 


17.352 
-15.642 
-3.952 


18.877 

-18.797 
-2.322 


20.161 

-21.486 
-0.917 


^^on— shelly 
^^off-shell^ 


-21.355 



-17.399 
1.757 


-18.797 



-19.407 

-2.079 



t without electromagnetic interactions 
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TABLE III. Deuteron Z)-state probability Pd, quadrupole moment Qd and asymptotic D/S 
ratio for various values of ji. 



fi Pd (%) Qd (fm2) D/S ratio 

(exact) 5.73 0.271 0.0258 

-1 (unitary) 5.16 0.268 0.0253 

-1 (exact) 4.98 0.266 0.0253 

+1 (unitary) 6.34 0.273 0.0263 

+1 (exact) 6.26 0.273 0.0260 
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^18 

Unitary trans, to 

^18 

Unitary trans, to \i=+l 

^18 



u/r 



w/r 



2 3 

r (fm) 



6.5 




Q Q Unitary transformation 



4.5 



-1 







1 



3000 




-150 I 

0.0 0.5 1.0 1.5 2.0 

r (fm) 



